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Recent studies show indication of the effectiveness of synchronization as a data assimilation tool for small or 
meso-scale forecast when less number of variables are observed frequently. Our main aim here is to understand 
the effects of changing observational frequency and observational noise on synchronization and prediction in 
a low dimensional chaotic system, namely the Chua circuit model. We perform identical twin experiments 
in order to study synchronization using discrete-in-time observations generated from independent model run 
and coupled unidirectionally to the model through x, y and z separately. We observe synchrony in a finite 
range of coupling constant when coupling the x and y variables of the Chua model but not when coupling 
the z variable. This range of coupling constant decreases with increasing levels of noise in the observations. 
The Chua system does not show synchrony when the time gap between observations is greater than about 
one-seventh of the Lyapunov time. Finally, we also note that prediction errors are much larger when noisy 
observations are used than when using observations without noise. 



The possibility of using chaotic synchronization 
for data assimilation is entertained. In particular, 
synchronization observed for two chaotic Chua 
models under unidirectional coupling is shown 
to be of relevance for forecasting using data as- 
similation techniques. Our results indicate that 
prediction method employed in the present is ex- 
tremely sensitive to observational frequency, Lya- 
punov exponent and the phase at which the cou- 
pling is switched off. 



I. INTRODUCTION 

Data assimilation (DA) is a powerful and versatile 
method for combining partial, noisy observational data 
of a system with its dynamical model, generally numer- 
ically implemented, to generate state estimates of non- 
linear, chaotic systems. A variety of data assimilation 
methods, broadly separated into deterministic or proba- 
bilistic methods, have been developed over the past few 
decades and used mainly in the earth sciences.'^El On the 
other hand, when two or more chaotic systems are cou- 
pled, they may adjust the properties of their motion due 
to coupling in such a way that their evolution becomes 
synchronized. Several types of coupling, unidirectional, 
mutual, common driving, etc., and their effects in the 
form of a variety of synchronizations such as identical, 
generalized, phase, lag, amplitude, etc., as well as appli- 
cations of synchronizations h ave been and continue to be 
a very active area of research.'^ 
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We will be interested in this paper with sequential 
data assimilation, which incorporates observations se- 
quentially in time. This method can be thought of as syn- 
chronization through unidirectional coupling, with obser- 
vations acting as the master syst em a nd the numerical 
model acting as the slave system.l^E^tl^l fact, this is 
exactly the nudging method for assimilation, suggested 
first in meteorology,^ - and later used for a number of dif- 
ferent studies.l^HIll Recent studies show that such tech- 
niques can also be used for assimilating data in order 
to track complex spatio-temporal dynamics of excitable 
media,!^ and forecast the state of a time-delayed high- 
dimensional system, e.g., in chaotic communication.!^ 
Synchronization based DA method may be useful for 
small or mesoscale predictions, specially when observa- 
tions are freq uent and coupled to a small number of 
variables .'i^I^ This method is relatively fast and robust, 
and can be implemented easily.'^il 

Two of the main characteristics of the observational 
data used in earth sciences applications is that they are 
discrete in time and are noisy. The main aim of this paper 
is to study synchronization with these characteristics in 
mind, as we explain below. 

1. It has been observed recentlj^^lHlIl that the time pe- 
riod between two observations is one of the impor- 
tant factors which significantly affect conventional 
DA schemes. Though frequent updates may pro- 
duce better results depending on the quality of the 
model as well as the observations, they increase the 
burden of computational cost.^Sl Since the choice of 
observational frequency depends on many factors 
such as cost and ease of observations, it is impor- 
tant to understand in detail how the time period 
between observations affects the accuracy of pre- 
dictions which in turn depends on the accuracy of 
the data assimilation. 
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For these reasons, we investigate the effect of vary- 
ing observational frequency, i.e., the gaps between 
the observations, on discrete time synchronization 
- this discussion forms the first part of the paper. 
As expected, we find that observations which are 
farther apart lead to poorer synchronization and 
consequently low accuracy of prediction. We also 
find that there is a certain threshold, which in the 
system we study is roughly of the order of the in- 
verse of Lyapunov exponent of the system, such 
that using observations with a time period greater 
than that threshold leads to no synchronization. 

It is of course well known that the observational 
noise plays crucial role in data assimilation as well 
as in synchronization. ^^^^^ Thus we also carry out 
the study described in the above paragraph with 
noisy observations with different noise levels, in or- 
der to investigate the effect of noise on synchro- 
nization. The discussion of these results forms the 
second part of the paper. Again as expected, the 
increase in noise leads to poorer synchronization 
and prediction. 



We present these results in detail in section [III] and a dis- 
cussion of these results along with directions for further 



studies are presented in section IV 



In order to investigate the effects of these two aspects, 
we perform the so-called "identical twin experiments"'^ 
as follows. We generate observations using a known tra- 
jectory, called the "true" trajectory of a numerical model. 
We use these observations (master) to try to synchronize 
the same numerical model starting with a different initial 
condition (slave). We quantify the "degree of synchro- 
nization" by calculating the root mean square (RMS) er- 
ror of the slave trajectory with respect to the true trajec- 
tory, leaving out the first half of the time-span, in order to 
get rid of transients. We assess the effect on predictions 
by calculating the RMS error over certain time-spans af- 
ter the synchronization is stopped and the slave model 
is let to run by itself. Currently we are investigating the 
effect of these two factors, the time period between obser- 
vations and the noise level, on data assimilation methods 
such as ensemble Kalman filter (ENKF).^ 

The dyna mical model we use is that of a chaotic Chua 
circuit .ISiEI] One of the main reason for this choice is that 
in future we would like to use the actual data from the 
circuit experiments being currently performed by two of 
the authors. •^^ This model is described in detail in the 
following section. 



II. CHUA CIRCUIT MODEL AND SYNCHRONIZATION 

To study synchronization we have used the Chua model 
used in Ref. [32l The model of the circuit is made dimen- 
sionless by substituting r = fit, where /? = (i?iCi)~^, 
and i?i and Ci are the resistance and capacitor used in 




FIG. 1. Two Chua attractors for the system in Eqs.|T|. In- 
set shows the large attractor, which represents the diverging 
dynamics, along with the above two attractors. 



the circuit respectively. The dimensionless model is given 
by. 



R 



dx _ 1200 

dy 
dr 

dz 



(y-x)~ g{x) , 



1200 

c[y + r'z] 



{y-x)+z 



(1) 



where 



-x \x\ < 1 

j{x) ^ { -[1 + b{\x\ - l)]sign{x) 1 < |a;| < 10 

10[(|a;| - 10) - (96 + l)]sign{x) 10 < \x\ 

and fj, c, r', and b are the parameters whose values de- 
pend on the capacitors, resistances and inductor used in 
the circuit .[23 The resistance R is the control parameter 
of the system. Depending on R, the system shows chaotic 
or regular behavior. In the present study, we have chosen 
a = 4.6/69 w 0.067, c = (1200 x 3300 x 4.6/8.5) x IQ-^ w 
0.779, / = 85/1200 w 0.071 and 6 = 1 - 1200/3300 w 
0.636. The model was integrated using the fourth or- 
der Runge-Kutta numerical scheme with integration step 
At — 0.01, and the output of the model has been ex- 
pressed in terms of real time t — t//3 rather than in 
terms of the dimensionless time r. 

The model equations show three attractors depending 
upon the initial conditions as shown in Fig.jljinset), and 
zoomed version of the smaller attractors by black and 
gray plots. As the basin of attractors is not known, it is 
not possible to tell a priory which initial condition will 
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lead to a particular attractor. Generally, for high val- 
ues of initial condition (approximately greater than 1 in 
absolute value), system goes to the large limit cycle at- 
tractor [inset of Fig. [l], and for smaller initial conditions, 
system stays on any of the two small attractors [Fig. [l] . 
As the larger attractor is not of experimental interest, 
the dynamics on this attractor is considered as divergent 
dynamics. 

In the present simulation, we have used R = 124517 for 
which the system shows chaotic behavior. The "true" (or 
master) initial condition is chosen such that the trajec- 
tory stays on any one of the small attractors depending 
upon the initial conditions. The largest Lyapunov expo- 
nent (A) is positive (« 4.04 bits/msec), and correspond- 
ing Lyapunov time is tx = 3.6 x 10~^sec. A was estimated 
using Rosenstein techniques-'^ and is consistent with the 
earlier results.!^ 

In order to generate any one set of observations, we 
first simulated a long trajectory on one of the smaller 
attractors. We refer to such a trajectory as "truth" in 
keeping with other data assimilation literature. Next, 
we sub-sampled this trajectory to extract several sets of 



discrete-time equispaced observations with varying time 
interval between the observations, which will henceforth 
be denoted by Tobs- For the present study, observations 
were generated with four different values of Tobs, namely, 
2.0 X 10"^ 5.0 X 10"^ 7.0 X 10"^ and 1.0 x 10"'' sec. 
Noisy observations were created by adding noise of par- 
ticular intensity (D) to these discrete-time observations, 
where D was chosen to be a specified fraction of the stan- 
dard deviation of the attractors. The results presented 
in Sec. |III| are for one specific set of observations, but 
the qualitative features of these results were found to be 
identical when we used observational sets generated from 
many different initial conditions. 



A. Master-slave coupling 

The model (slave) was unidirectionally coupled with N 
discrete-time observations (master) at times ri , . . . , tjv as 
follows. 



dx 1200 



N 



dr R 



{y- x) - g{x) - kx'^ [a;(r) - Xobsin)] S{t - n 



dy 
dr 

dz 
dr 



1200 



{y-x) + z 



1=1 

N 



N 



-c[y + r'z] - fc^ ^ [z{t) - Zobsin)] 5{t - n 



(2) 



where, kx,ky, and kz represent the x, y, and z coupling. 
In the case of continuous time observations which is also 
discussed below, the slave is coupled to the master di- 
rectly, by replacing the summations of delta functions in 
the above equations by terms such as kx[x{T) — a;obs(''')]- 
We have quantified synchronization using RMS error of 
the difference of the synchronized and true trajectory 
over the last half of the trajectory, and the quality of 
prediction is again quantified by the RMS error with re- 
spect to truth over varying time periods, as discussed 
in detail below. The results about synchronization and 
prediction for different Tobs and noise level (D), and a 
comparison of these results with the case of continuous 
time synchronization are presented in the next section. 



III. RESULTS 

The main aim of this section is to discuss the effects 
of the observational time period Tobs and the observa- 
tional noise variance D on the root mean square of the 
difference between the slave system and the master, both 



during the time period in which the slave is coupled to the 
master (synchronization phase) and after the coupling is 
switched off (the prediction phase). 



A. Synchronization with observations without noise 

(D = 0) 

In order to investigate the performance of the synchro- 
nization, we have first studied the case of coupling only 
the a:;-component, i.e. setting ky = = kz in Eqns. [2] 
while kx is varied from zero to higher values. We quan- 
tify synchronization of each components by calculating 
the RMS error of the component-wise difference of mas- 
ter and slave trajectories, and also the total RMS er- 
ror. For continuous case, the slave model synchronizes 
with the master within a very short time. Solid lines in 
Fig [2] (a) and (b) show the RMS errors in x-component 
and total errors respectively. They show that the system 
started synchronizing around kx ~ 0.005, and achieved 
full synchronization around kx ~ 0.03 as the RMS errors 
are almost zero. The system remained in synchronized 
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FIG. 2. RMS errors in the a;-component (top) and the total 
RMS error (bottom) when a;-component was observed. In 
both these plots as well as in Fig. [4j thick solid line is for 
continuous time synchronization while dashed, thin solid, and 
gray lines are for cases with Tobs = 2.0 x 10~^, 5.0 x 10~^, 
and 7.0 x 10~^ sec respectively. The y- and ^-components 
also show similar qualitative behavior. Note that RMS error 
above approximately 0.5 indicates loss of synchronization. 



state even at higher values of k,j. . 

Synchronization with x-component observations at dis- 
crete time is shown in Fig. |2] In this figure, the dashed- 
lines show the RMS errors when ^-observation were 
recorded with Tobs = 2.0 x 10~^ sec, which is approx- 
imately 18 times less than t\. The figure also shows 
that the model started synchronizing around — 0.18 
with the observations and achieved full synchronization 
around > 0.5. For k^ > 2.5 model diverges from the 
observations. So the range of the coupfing k^ for which 
there is synchrony has become finite when observations 
were taken at a finite interval. Synchrony is also clear 
from the plots of x and y time series of the observational 
trajectories and slave trajectories shown in Fig. [3j They 
show the model (gray lines) gets perfectly synchronized 
with the observational trajectory (black lines). 

Note that RMS error above approx. 0.5 indicates a 
loss of synchronization, hence effectively synchronization 
is obtained only for the first two of these cases, i.e., con- 
tinuous time observations and with Tobs = 2.0 x 10~^ 



sec. 



Black thin solid lines [Fig [2] show the RMS errors, 
when the observations were taken with Tobs = 5.0 x 10~^ 
sec. {fa ^t\). In this case errors decreased slowly with fcj: 
and showed partial synchronization only in a very small 
range of k^ around k^ « 1.7. The time series at fc^; — 1.74 
in plots (c) and (d) of Fig. [s] show such partial synchrony 
between the observed trajectory (gray lines) and the slave 



FIG. 3. First and second columns represent the x and y com- 
ponents of the slave trajectories (solid line) and model trajec- 
tories (gray line) at k^^ = 1.74 for the observations taken at 
Tobs = 2.0 X 10"^ (top) and 5.0 x 10"^ (bottom). Note that 
the jumps in a;-component and in j/-derivative are because of 
the (5-functions in Eqs. ([2|. 



model (black line) for Tobs = 5.0 x 10~^ sec. The be- 
haviour of the z-component was similar to that of the y. 
For Tobs = 7.0 x 10~^ sec, the RMS error is least around 
kx = 1.19 but even in this case, it is far from being syn- 
chronized. These results show that Tobs is an important 
factor in synchronization that must be considered in DA 
experiments. 

When y-component was observed, i.e., when we set 
kx — Q ~ kz and vary ky, the model showed synchrony 
for different range of ky than that was for x observa- 
tions. In Fig. |4j the solid lines show the RMS errors in 
(a) a;-component and (b) total RMS error for continuous 



observation case. Around fc,. 



0.2 the model achieved 



synchrony and remains in this state even at higher val- 
ues of ky. For Tobs = 2.0 x 10~^ sec (dashed lines), slave 
system achieved full synchronization around ky = 0.53, 
and remains synchronized for 0.53 < ky < 2.22. When 
Tobs = 5.0 X 10~^sec (thin solid), the synchronization re- 
gion shifts to 1.14 < ky < 2.72. Similar to x-observation 
case, there is no synchronization for the y-observation 
with Tobs — 7.0 X 10^^ sec (gray line) and above. Thus 
we see that the y-component is a more dominant dynam- 
ical variable than the x-component. 

When model was coupled with master through z- 
observations {k^ — — ky), there was no synchronization 
even for the continuous observations case. Panel (a) of 
Fig ^ shows the RMS errors in x (solid line), y (dashed 
line)7and z (dotted line) components as kz is varied. This 
is also clear from the plots of cc, y, and z time series of 
the master (gray) and slave (solid) shown in panels (b). 
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FIG. 4. RMS errors in the x-component (top) and the total 
RMS error (bottom) when a;-component was observed. The 
lines are for same Tobs as described in Fig. [2] Comparing 
the RMS error for Tobs = 5.0 x 10"^ sec (thin sohd lines) in 
this figure and in Fig. [2] we see that y-observations are more 
informative than the ^-observations. 



2 



0.5 




(a) 



0.5 1 1.5 2 2.5 

coupling constant 




FIG. 5. (a) RMS errors in x (solid), y (dashed) and z- 
component (dotted line) when z-component was observed 
continuously in time, showing that there is no synchronization 
when z is observed. Time series of slave (gray) and master 
(black line) at = 1.2 is shown in (b) x, (c) y and and (d) z. 
The X- and y-components of the slave have been multiplied 
by factors of 20 and 0.5 respectively after subtracting mean, 
in order to see the comparison clearly. 




FIG. 6. Attractors of the master (solid black) and slave sys- 
tem (dotted line) when z was observed. The structures of 
these attractors are clearly different from each other. 



(c), and (d), respectively. The difference between the 
master (solid line) and slave (dotted line) is also clear 
from the attractors shown in Fig |6] A possible expla- 
nation for not having synchronization with z coupling 
is that this variable may not contain sufficient informa- 
tion about the dynamics. Hence in the current regime 
of the Chua system x and y contain useful information 
for synchronization but not the z-component.li^ Finally 
as the model did not show any synchrony for continu- 
ous case, obviously we did not get any synchrony with 
discrete time observations either. 



B. Synchronization with noisy observations 

To study the effect of noise on synchronization, we 
added noise to the x-component of observ ation with 

-^obs 
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and 



2.0 X 10~ sec, which is used in Sec. 
repeated the same experiments. We have not performed 
any experiments with noisy continuous time observa- 
tions, since we are only considering deterministic models. 

In Fig.[7](a) and (b), the solid lines represent the RMS 
errors in x-component and total error respectively when 
added noise intensity \%D — 2%. In this case, we got syn- 
chrony in the range of 0.66 < < 2.1, and the model 
started desynchronizing around = 2.1 as shown in 
Figs.[7] So the noisy observations squeeze the k range fur- 
ther. With increase in D overall range remained almost 
same. In Figs. [7| solid, dashed, dotted, dashed-dotted, 
and gray lines are the RMS errors when D = 10%, 20%, 
30% and 40% respectively. 

Effect of the noise is also clear from the time series 
plots. Fig. [8] shows the time series of the observation 
(thin solid line) and slave trajectories for D = 10% (thick 
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FIG. 7. RMS error of synchronization when noisy x- 
component observations with Tobs ~ 2.0 x 10~^ are cou- 
pled to the slave model: (a) RMS error in x-component and 
(b) total RMS error. The solid, dashed, dotted, dashed- 
dotted, and gray lines represent the RMS errors for noise 
levels D = 2%, 10%, 20%, 30%, and 40%, respectively. 
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FIG. 9. Prediction RMS errors in the a;-component for dif- 
ferent tfh when X was observed without noise with Tobs = 
2.0 X 10"^ (top) and Tobs = 5.0 x lO"'^ (bottom). Plots in 
the right column are the zoomed version of the boxed regions 
of the left plots. In each plot, the thick solid, gray, dotted, 
dash-dot, and thin solid lines show the prediction RMS errors 
when tfh — jQ^x, t\, 5tA, Wt\ and 50t\ respectively. 
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time (msec) 

FIG. 8. Time series of the observed trajectory (thin black), 
and the slave trajectories coupled to observations with D = 
10% (thick gray) and 40% (thick black) noise with Tobs = 
2.0 X 10~^sec. The noisy observations with D = 40% are also 
shown (stars). 



gray) and 40% (thick solid line). They show that the 
deviation of the slave trajectories from the observed one 
increases with increasing noise level. But z-component 
deviate lesser than x and y-component (not shown). It 
is also observed that the RMS errors increases almost 
linearly (not shown here) with increases in D. 



C. Prediction with truths and noisy observation 

In this section we consider the use of the slave model 
for prediction. This is done in the following manner. We 
couple the observations with the slave system through 
Eqs. ([2]) for a long time period, in particular, until tcoup = 
10 msec in our numerical experiments. Then we switch 
off the coupHng, i.e., use Eqs. ([I]) for the slave beyond this 
time. The RMS errors described below are calculated 
over varying time periods, which are called the forecast 
horizons tf^- We present the results for x-observations 
with Tobs = 2.0 X 10"^ and Tobs = 5.0 x 10~^ since larger 
Tobs do not show any synchronization. 

In Fig. [9] the prediction RMS errors in x-component 
has been shown for the range oi < < 5. The y- and 
z- components follow the same qualitative trend. Thick- 
solid, gray, dotted, dash-dotted and thin solid lines in 
Fig. are the RMS errors for tfh = jfjtx,tx,5tx,l0tx, 
and bOtx respectively, where boxed regions of Fig. [9] 
(a) has been zoomed in Fig. |9] (a'). We see that for 
tfh = jQtx,t\,5t\,i0tx slave system predicts well, i.e.. 
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FIG. 10. Prediction RMS errors in the x-component for differ- 
ent tfh when X was observed with noise levels D = 6% (top) 
and D = 10% (bottom) with Tobs = 2.0 x 10"^ Plots in 
the right column are the zoomed version of the boxed regions 
of the left plots. In each plot, the thick solid, gray, dotted, 
dash-dot, and thin solid Hnes show the prediction RMS errors 
when tfh — jxi^x, tx, 5t\, Wt\ and 50tA respectively. 
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FIG. 11. Differnce between the observed and the predicted 
trajectories, when the coupling is switched off at f = 10 msec. 
Top panel (a) is for the case of observations with no noise, 
and the bottom panels (b) and (c) are for the case with noisy 
observations with intensity D = 6% and D = 10% respec- 
tively. 



the RMS errors are small, but when tfh > I0t\ predic- 
tion diverges from true trajectory. Figs [o] (a') [thin solid 
line] show that RMS error at tj^ = 50tx is significantly 
large. Fig. 11 (a) show that slave model becomes com- 
pletely uncorrelated with the master around tfh — 40<a- 

Another important observations from our numerical 
experiments is that the forecast horizon for which the 
slave shows small RMS error with respect to the master 
is highly dependent on the state of the system at which 
the coupling is turned off. This is because the uncer- 
tainty growth rate is different on different parts of the 
attractor.l^ 

When noisy observations are coupled with the slave 
system, the divergence of the slave from the master is 



rapid once the coupling is switched off. Fig. 10 shows er- 
rors in cc-components when noise with intensity D — 6% 
(top panels) and 10% (bottom panels) is added to the ob- 
servations. We have seen that for tfh < \tx (not shown 
in figure), there is a range of coupling constants for which 
the prediction RMS error is small enough to give good 
prediction. For larger tfh^ the slave model become un- 
correlated with the master and predictions cannot be ob- 
tained. This is the case for both cases of the observations 
with 6% and 10% noise. Divergence is also clear from 
the difference of ccobs ~ x time series of with 6% of noise 



[Fig. 11 (b)] and 10% of noise [Fig. nil (c)] respectively. 



IV. CONCLUSION 

This paper is devoted to the study of discrete-time 
unidirectional synchronization, with specific emphasis on 
aspects which are of relevance to data assimilation prob- 
lems. In particular we consider the observations of a 
system, in this case the chaotic Chua system of Eqs. ([I]), 
to be the master coupled unidirectionally to the slave 
system of Eqs. ([2| and examine how the slave system 
synchronizes with the master as we change various char- 
acteristics of the observations. The main discussion is 
focused on studying the effects of changing observational 
period, i.e. the time between the observations Tobs: £^nd 
observational noise levels D. Some of the main results as 
well as directions for future research which we are cur- 
rently pursuing are discussed below. 

We obtained synchronization with continuous time ob- 
servations of the X- and y-components of the system, but 
not with 2-component observations. When discrete ob- 
servations with finite Tobs were used, Chua model showed 
synchrony for a finite range of coupling constants, and 
this range decreases with increase in Tobs and noisy ob- 
servations. The study also shows that the system shows 
synchrony only when the observations were taken with 
Tobs much less than t\ (almost ^t\) while for larger Tobs 
there is no synchronization between the master and the 
slave. 

In many cases of practical interest in earth sciences. 
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data assimilation is used as a tool to improve prediction 
of the observed as well as unobserved components of the 
system. With this in mind, we also study the predic- 
tion errors over different forecast horizons after the cou- 
pling between the observations and the slave is turned 
off. We see that in the cases when synchronization errors 
are small, the prediction errors even for forecast hori- 
zons far beyond the Lyapunov timescale are small. This 
happens only for the case of observations which are not 
noisy. When using observations even with a very small 
noise, the slightly higher synchronization errors coupled 
with chaotic nature of the system lead to loss of predic- 
tive power over forecast horizons which are a fraction of 
the Lyapunov timescale. These results indicates that the 
association of the Lyapunov exponent with predictability 
horizons has limited use."^^ These results also show the 
limitations of discrete-time synchronization with noisy 
observations and the importance of the study of other 
data assimilation techniques which systematically take 
into account the noise in the observations. 

In synchronization choice of coupling constant k is 
merely ad hoc and depends on a particular application. 
In case of data assimilation using similar techniques, e.g. 
nudgin^^, the choice of k can be interpreted in terms of 
variational approach. In particular, Ref. [18] shows, us- 
ing the variational approach, that the optimal coupling 
constant is inversely proportional to observational noise 
covariance: k ^ R^^, where R is the covariance matrix 
of observational errors. As the range of k on which model 
shows synchrony depends mainly on observed component 
and noise, we are investigating the possibility that a sim- 
ilar formulation may also be possible for synchronization 
to choose the appropriate k. We will also be comparing 
the synchronization based methods with data assimila- 
tion techniques such as the ensemble Kalman filter^ and 
Bayesian methods.'^ 

Another ongoing investigation is to study the perfor- 
mance of these methods in the chaotic Chua model using 
real data from the circuit!^. This study will give us an 
insight about the choice of coupling constants and obser- 
vational time period, the effects of observational noise as 
well as of the errors in the numerical model. 
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